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We derive and solve flow equations for a general 0(A^)-symmetric 
effective potential including wavefunction renormalization corrections 
combined with a heat-kernel regularization. We investigate the model 
r~| , at finite temperature and study the nature of the phase transition 

in detail. Beta functions, fixed points and critical exponents f3, v, 5 
and r] for various are independently calculated which allow for a 
verification of universal scaling relations. 



1 Introduction 

The manifestation of universality of critical phenomena enables the applica- 
bility of the massless 0{N) scalar field theories to a wide class of very differ- 
ent physical systems in the vicinity of the critical temperature Tg. E.g. for 
= 4 this theory is used as an effective model (the linear a-model) for the 
chiral phase transition in two light-quark flavor QCD which may have the 
same universality class [1]. On the other hand in condensed matter physics 
the case = 3 corresponds to the well-known Heisenberg model describing 
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the ferromagnetic phase transition. Further apphcations (hke the Kosterhtz- 
Thouless phase transition (A^ = 2), hquid- vapor transition describable by 
the continuous spin Ising model {N = 1) or even the statistical properties of 
long polymer chains (A^ = 0)) are often discussed in the literature. 

In the vicinity of a second-order phase transition the equation of state 
(EoS) obeys a universal scaling form, which cannot be assessed by ordinary 
perturbation theory due to infrared divergences. For a negative mass term 
the A^-component theory exhibits spontaneous symmetry breaking of the 
0(A'')-group down to 0{N — 1) which is restored at sufficiently high temper- 
ature. According to the Goldstone-theorem {N — 1) fields (the pions) become 
massless in the chiral limit and the remaining one (the radial sigma mode) 
stays massive. Due to these massless modes a naive perturbation expansion 
for the effective potential fails to give an appropriate description and the 
potential itself develops an imaginary part around the origin. By means of a 
nonperturbative Wilsonian renormalization group approach one can circum- 
vent these difficulties and go beyond any finite order in perturbation theory. 

In this work we present a renormalization group (RG) approach formu- 
lated in the Euclidean space. It is combined with a heat-kernel regularization 
which allows a straightforward calculation of the effective potential with an 
0(iV)-symmetry for arbitrary N, not only at finite temperature but also di- 
rectly at the critical temperature. We determine several critical exponents 
and investigate the phase transition regime in detail going beyond the so- 
called 'local potential approximation' to also include wavefunction renormal- 
ization corrections. 

The paper is organized as follows: in Sect. 2 we summarize the Wilsonian 
renormalization group concept and show how fiow equations combined with 
a proper-time (heat-kernel) regularization are obtained for a general 0{N)- 
symmetric potential. Based on a perturbative one-loop effective potential 
expression we derive flow equations for the potential and wavefunction renor- 
mahzation for a d-dimensional massless 0{N) scalar theory. These coupled 
RG flow equations are further improved by taking into account the contin- 
uous feedbacks from the higher modes to the lower ones. The derivation of 
the flow equations with wavefunction renormalization corrections is given in 
Sect. 3. 

Sect. 4 details the numerical results of the solution of the flow equations. 
We start with a discussion of the numerical implementation of the properly 
rescaled equations for a general potential with 0(A'^)-symmetry on a grid 
whereby the full potential is numerically taken into account without resorting 
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to polynomial approximations. The order parameter and the potential itself 
are studied at finite temperature. 

Sect. 5 is dedicated to the critical behavior of the system at the transi- 
tion temperature. The scaling solution, fixed points and beta functions are 
investigated in detail and numerical calculations of the critical exponents for 
different N including the anomalous dimensions are presented. They are 
found to be in perfect agreement with other works and approaches. Sect. 6 
contains the summary and concludes our discussion. Technical details are 
relegated to the appendix in order to improve the readability. 

2 Concept of the RG method 

The effective action r[$] is defined by a Legendre transformation of the gen- 
erating functional VF[J], which generates all connected Feynman diagrams 
and can be expressed by a functional path integral. Except for numerical 
simulations, these integrals usually cannot be calculated exactly. In addition 
they are plagued by infrared (IR) divergences. Instead of performing the 
integration in one step one can follow the idea by Wilson and Kadanoff [2] 
and split the field $ into fast-fiuctuating and slowly-varying modes 



and perform the functional integration only over the fast-fiuctuating compo- 
nents with momenta q larger than the scale k separating both field modes 
from each other. The irrelevant short-distance modes which are insensitive 
to the physics at large correlation scales are thus ehminated. One is left with 
a low-energy "effective action" parameterized by the averaged fields at the 
IR scale k which acts as a momentum cutoff. The full effective action de- 
scribing the one-particle-irreducible (IPI) propagators and vertex functions 
is generated by the renormahzed effective action in the limit A; — > 0. Thus 
the /c-dependent action provides a smooth interpolation between the bare 
(classical) action at an arbitrary ultraviolet (UV) scale A where no quantum 
ffuctuations are taken into account and the renormahzed action in the IR. 
For further details we refer the interested reader to refs. [3], [4], [5], [6]. 

The effective action is characterized by infinitely many couplings mul- 
tiplying all possible invariants which are consistent with the symmetries of 
theory in question. In order to proceed one has to make approximations. 
One possibility is the expansion of the effective action in powers of the fields 




(1) 



3 



about an arbitrary field ^o{x). The Taylor coefficients of this expansion 
r*^"'^(a;i, .... Xn) then represent the full IPl propagators and vertices. 

Alternatively one can expand the effective action in powers of derivatives 
corresponding to an expansion in powers of the momenta which we will pursue 
in this work.^ The invariants are in this case classified by the number of the 
appearing derivatives and the coefficients are functions of constant fields. 
For a one-component scalar theory in d-dimensional coordinate space this 
expansion takes the form 

r[$,] = /d'^x{-n<^'c) + ^^(<^'9(W' + m')((W')' + --j • (2) 

The lowest-order term, the effective potential V, should be a convex function, 
with no dependence on derivatives and can be rewritten as a sum of all 
IPl Green functions r*^"-' with n external lines carrying zero momenta. In 
other words the effective potential is the generating functional of the zero- 
momentum IPI Green functions. In the following the expansion (2) will be 
truncated at second order (9(9^) to derive a set of coupled partial differential 
equations for the effective potential V and wavefunction renormalization Z 
at the scale k. 

3 Application to the 0(A^)-model 

In this section the RG method will be apphed to the 0(A^)-symmetric linear 
sigma model. We concentrate on the meson fields (A^ — 1 massless pions and 
one massive radial sigma mode) and work with an arbitrary N instead of 

= 4 fixed. 

We employ the Euclidean metric due to the heat-kernel regularization for 
the real part of the functional determinant^. This metric also allows for a 
straightforward extension to finite temperature within the imaginary time 
formalism as discussed in Section 3.3 ^. 

^Which expansion pattern for the effective action should be used depends on the 
underlying theory. The derivative expansion will not be useful in situations where the 
momentum-dependence of the correlation functions are manifestly influenced by the in- 
teraction. 

^For notations and details see e.g. [7] 

^An implementation of the Wilsonian RG in the real-time formulation can by found 
e.g. in [8] 
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At the scale A in the uhraviolet region of the theory we define the effective 
Lagrangian for the 0(A^)-model by 



U = I {d,$y + V{$') ; = ^ - ^if (3) 

with the A^-component vector $ = [^1,^2, ■ ■ ■ ,^n)- The negative sign in 
the potential V signals the broken phase and a positive sign indicates the 
symmetric phase. 

The field $0 denotes the minimum of the potential. If this field is con- 
stant only the lowest-order term (the effective potential) of the derivative 
expansion of the effective action Eq. (2) can be computed. In the literature 
this approximation is called the 'local potential approximation'. 

Higher-order terms of the derivative expansion can be calculated with a 
varying x-dependent background field ^o{x) (see e.g. [9], [10], [11]). The 
expansion of the effective action of the 0(A^)-symmetric model up to order 
C(a^) reads 

m = J d'x {-l^($2) + ($2) (5^$)' + ($') . (4) 

To this order there are two independent wavefunction renormalization terms 
Zi and Z2 for AT > 1 

In the following we will sketch the derivation of the flow equation in- 
cluding the wavefunction renormalization corrections. The perturbative one- 
loop contribution to the effective action yields formally a non-local logarithm 
which must be regularized. Proper-time regularization is the most commonly 
used regularization for the real part of the determinant and results in a finite 
local action (see e.g. [12]). It was introduced by Schwinger and maintains a 
finite sharp dimensionful proper-time cutoff [13]. We replace this sharp ultra- 
violet cutoff for the proper-time variable r by a smooth, a priori unknown, 
kernel or 'blocking function' fk in the integrand which will be specified later. 

After introducing a complete set of plane wave states in the heat kernel 
[14] this results in the following d-dimensional effective action 

r[$] = -\j d'x f J ^tr e-Kp-^wa.-a^.-^:;w) (5) 

*The next order in the derivative expansion would yield ten independent terms. 
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with the shorthand notation V-j{^) — for the {N x A'")-matrix valued 
second derivative of the 0(A'")-symmetric potential V given by 

y/' = A - ^l) 6,, + 2A$,$,- . (6) 

The trace in Eq. (5) runs over the N fields and can be computed analytically 
with standard techniques [15], [16]. In the following we omit the index struc- 
ture of the potential term V^j and follow the techniques outlined in ref. [17] 
which we generalize to an 0(iV)-symmetric matrix- valued potential V". 

The expansion of exp [— r {—2ip^dn — d"^ -\- V")] in Eq. (5) in powers of 
derivatives up to second order yields 



{oo /T///\™ n oo f 1 \n— 1 n n— 1 

n=0 n=2 fe=0 

oo / 1 \n— 1 n n— 2 n—2—k ~j 



(27r)'^ 

oo / 1 \ra — 1 _ra n — 1 

n— 1— fc 



n=3 k=0 1=0 ) 

where all linear terms in the momenta are ignored since they vanish any- 
how after odd p-integration. The above expression can be further simplified 
by a tedious but straightforward calculation ^. All infinite sums converge and 
can be evaluated analytically with careful attention of the non-commuting 
operator order. 

To improve readability we will use in the following the abbreviations 

= A(3$2 - <1)^) and = A($2 - (8) 

and exploit the 0{d) symmetry of the momentum integration. 
To order 0{d*) the effective action splits into two parts: 

with r*^*') denoting the one- loop potential contribution (no derivatives) 

r'») = -1 / A / 1/<»' / ^e--' {e-< + (N - l)e-'<} (9) 
^For technical details see the appendix A. 
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and r^^^ the second-order contribution (containing two derivatives) 



{2n) 




d,..mld„rnl e"^'"' 



+ {N-1) i^d'ml - ^d,mld,ml] e 



ml — ml 



2 3 



($2)5 



$2 



2t^ T-4 \ 

[—d'^ml - —d^mld^mlj e 

+(N - 1) [^av 



(10) 



+2t 



( — ^ 



$2 



($2)2 



For iV = 1 this outcome for F'^^^ agrees with that obtained in ref. [17] using 
a similar procedure. 

In order to extract the two different wavefunction renormahzation Zi 
and Z2 of the effective action F in Eq. (4) one has to collect the terms in 
front of (dfj,^)^ and of ($9^$)^. To simplify the amount of work drastically 
we will use the truncation to a uniform wavefunction renormahzation which 
means wc neglect the field and momentum dependence of the wavefunction 
renormahzation. This approximation corresponds to considering only the 
wavefunction renormahzation constants at the minimum of the potential $0 
and neglecting their derivatives. 

For iV > 1 one has to work with two different Zi and Z2 correspond- 
ing to two anomalous dimensions. These are related to the wavefunction 
renormahzation for the Goldstone modes Z^^ and the radial mode Z„ via 



Ztt — Z\ 



As will be seen later the difference between the corresponding anomalous 
dimensions at the scaling solution vanishes. We can further simplify the 
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calculation by neglecting Z^- This approximation describes all the qualitative 

features at the phase transition without loss of quantitative predictive power. 

Comparing the expressions in (9) and (10) with Eq. (4) we can now extract 
the effective potential contribution 

and the Zx = Z wavefunction renormalization contribution 
1 Z c dT 1 



Z = 



/?7^[(/r-2/r) (12) 




where momenta have been rescaled by the corresponding wavefunction renor- 
malization. 

Note, that we have introduced two different blocking functions, jf^ and 
to be discussed in the next Section. To proceed we differentiate Eqs. 
(11), (12) with respect to the scale k which is introduced via the blocking 
functions thus obtaining the desired flow equations for the potential and 
wavefunction renormalization. In doing so one has to know the derivative of 
the blocking functions j'^ i = 0, 1, . . . with respect to k. 



3.1 The regulating blocking functions 

In this Section we sketch the way how to implement the regulating blocking 
functions f^ij) in the proper-time integrations^. They serve to cut off the 
diverging momentum contributions to the flow equations (see e.g. [17, 18]). 
Going beyond the local potential approximation by including higher contri- 
butions to the effective potential through the wavefunction renormalization 
we have to deal with different powers of momenta. In deriving flow equations 
for an ©(A'^) -symmetric potential we there encounter additional momentum 
integrals which force us to modify the blocking functions. In fact, at each 
level of the gradient- or derivative expansion new blocking functions must be 
taken into account (cf. e.g. [18]). In order to find a condition on the structure 
of the blocking functions we compare the operator regularization method to 

^Other authors [18] call the kernel /fc(r) also smearing function. 
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a momentum regularization method. In ref. [4] a differential equation for the 
blocking function in the local potential approximation was found by compar- 
ing the operator cutoff to a sharp momentum cutoff. This yields RG flow 
equations with a characteristic logarithmic structure similar to the Wegner- 
Houghton equation [19]. When the momenta are rescaled as — > Zp^ we 
obtain a similar d-dimensional equation^ 

The function Z^^-* regulates momentum integrals of the form 

J^e-^'^' (14) 
(27r)<i ■ ^ ' 

and plays a similar role as the covariant Pauli-Villars cutoff. Its derivation 
can be found in detail in refs. [4], [5]. 

The above expression must be modified when higher derivatives are taken 
into account because one encounters higher powers in the momentum inte- 
grals of the form e.g. 



Comparing the sharp momentum cutoff version of I2 with the proper-time 
expression with a new blocking function f^'^ and by taking the derivative 
with respect to k one obtains the following equation 

Thus each power of in the numerator of the integrand (cf. Eq.(15)) gener- 
ates an additional power in the proper-time variable r thereby accelerating 
the convergence of the flow equation. This result is identical to that in 
ref. [18]. 

For completeness we quote the result for the blocking function which 
regularizes momentum integrals of the form 



(27r)« 



e ''^'^^P,,PvP\Pp ■ (17) 



prime on /^'•' implies differentiation with respect to the scale k. 
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Although they are not needed in our approximation we will employ such 
higher blocking functions when we analyze the influence of the blocking func- 
tions on the convergence of the flow equations in Sect. 5.4. 

Comparing in the same manner with a sharp cutoff version yields the 
differential equation 



(2) / 



' d{d + 2)V{d/2) 



(18) 



3.2 The flow equations 

With the general differential equation for the blocking function f^^ at hand 
we can immediately obtain the flow equations. It turns out that the lowest 
blocking function f^"^ is sufficient to regulate the potential equation alone^ 
(no derivatives) but not for higher momentum contributions. The conver- 
gence of the solution of the flow equation is very poor with /|°^ when wave- 
function contributions are taken into account. This result agrees with that 
found in [4], [5]. Thus in order to accelerate the convergence from the very 
beginning we chose a blocking function f'^^ with i > 0. We start with i — 1 
and replace in Eqs. (11), (12) f'f' by Performing now the convergent 

proper-time integration this yields the flow equations 



dfV = S. 



for the potential and 



'It 



+ 



N -I 



1 + 2^;' + A^'^v" 1 + 2v' 



(19) 



Z 



2Sd 
^^Zk^ d 



+ 



1 + 4$2^//-)2 2$2-y/' V 1 + 4$2-y" 



1 



(20) 



for the wavefunction renormalization where we have introduced the notations: 
dV , ^ 2 



dtV = k- 



dk 



Zk^ 



and 



r((i/2)(47r)'^/2 



A prime on the potential V denotes differentiation with respect to $2. Due 
to our uniform wavefunction renormalization approximation (no derivatives 
of Z) we have to evaluate the flow equation at the potential minimum $o- 

*This blocking function corresponds to the cutoff function of type (II) in ref. [4]. 



10 



The functions in squared brackets describe the threshold behavior of massive 
modes and are sometimes called 'threshold functions' [6]. Each massive mode 
decouples from the evolution towards the infrared limit /c — > (see details 
in [5]). One advantage of our choice of the blocking functions is that the 
threshold functions can be obtained analytically. In Eq. (19) one recognizes 
the persistent contribution of the {N — 1) massless Goldstone particles to the 
evolution since v' vanishes for the minimum of the potential (cf also Eq. (20)). 

In the above flow equations we, in addition, perform a renormalization 
group improvement thus going beyond the independent-mode approxima- 
tion. The improvement consists of replacing the bare potential V (and its 
derivatives) on the rhs of the flow equations with the running A;-dependent 
counterpart V^. Then higher graphs such as daisy and super-daisy diagrams 
are taken into account. Thus we substitute 

ml = 2V^ + A^'V^' and ml = 2Vl , (21) 

which should not be conftised with the abbreviation in Eq. (8). Of course, 
the potential Vk is still a function of the field <I>^ itself. This type of nonlinear 
flow equations has to be solved numerically on a grid as discussed in the next 
Section. 

3.3 Finite temperature 

Here we wish to discuss the generalization to finite temperature. Within the 
Matsubara formalism this is easily accomplished by replacing the d-dimen- 
sional momentum integration in Eqs. (9) and (10) by a (li — l)-dimensional 
integration and by introducing a discrete Matsubara summation with bosonic 
frequencies a;„ = 2n7rT in the zero- momentum component. An additional 
length scale (the inverse temperature) is now introduced and besides quan- 
tum fiuctuations additional thermal fiuctuations appear. At finite temper- 
ature the scale k serves a generalized IR cutoff for a combination of the 
three-dimensional momenta and discrete frequencies. The influence of the 
different choices of blocking functions on the two types of modes can be 
found in more detail in [20]. We use in this work cutoff blocking functions 
with a (i- dimensional momentum variable. 

In order to be consistent with the zero-temperature flow equations we 
make the same approximations and use the same blocking function fj^^^ at 
flnite temperature. The finite-temperature extension of Eqs. (19) and (20) 
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with a uniform wavefunction renormalization contribution is 



dfV 



8(47r) 



MO, 



ZA;2' 2 



+ 



{N - 1)M ( 0, 



\ ' ZF' 2^ 



(22) 



where the foUowing abbreviations for the finite-temperature threshold func- 
tions 



2\p 



(23) 



have been introduced together with the short-hand notation of Eq. (21). 
For the uniform finite-temperature wavefunction flow equation we obtain 
accordingly 



dtZ 



Z Te (5 r(3/2) 



(47r)3/2 [2 ml -ml 
5r(5/2) 
4 

r(5/2) 1 



+ 



M 0, 



ml — ml 

r(7/2) 1 



-M l, 



2Zk^ A;2 



-M l, 



m 



Zk^' 2 y 
ml 5 

2 7\ 



+ M{0, 



Zk^' 2 

2 

mz 



.2 5^ 



-Ml, 



ZA;2' 2 



2 5^ 



ZF' 2 



+ 1, 



Zk^' 2^ 
ZA;2 ' 2 , 



mz 



(24) 



It is possible to calculate the zero-temperature limit of the finite-tempe- 
rature flow equation (22) and (24) analytically by using the zero-temperature 
limit of the threshold functions^ : 



M (0,m^,a) 



k 2 ■4- ■■(2a -3) 



1 



and 



M (l,m^,aj 



ttT 1 • 3 • • • (2a - 2) (1 + m2)(2«-i)/2 

A;^ 2-4---(2q;-5) 1 1 

27rT 1 • 3 • • • (2q; - 4) (a - 1) (1 + m2)(2°-3)/2 " 



(25) 



(26) 



In the next section we discuss the numerical solution of these equations 
for the broken phase. 



^see e.g. ref.[5]. 
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4 Results 



4.1 Numerical implementation 

We discretize the field $^ for a general potential term on a grid and 

numerically solve the coupled set of flow equations (19), (20) for zero temper- 
ature and (22)-(24) for finite temperature at each grid point. 

The flow equations incorporate derivatives of V with respect to the field 
up to second order. Thus further differentiation of the flow equation with 
respect to the fields increases the order of derivatives of V. Following the 
method described in [21] we calculate the first- and second derivative of the 
flow equation for the potential with respect to $^ and thus generate deriva- 
tives of V up to fourth order. The unknown derivatives V" and V"" are 
determined by matching conditions at intermediate grid points in order to 
close the highly coupled system of differential equations. This matching con- 
dition is based on a Taylor-expansion of V and V" up to forth order at each 
grid point $^ and it requires continuity of these expansions at intermediate 
adjoining grid points. 

For n grid points one obtains 2n — 2 equations for the 2n unknown deriva- 
tives V" and V"". The missing two conditions of the boundary grid points 
are determined by an expansion of the third derivatives. This results in an 
2n-dimensional algebraic coupled set of flow equations for V and V" which 
can be solved by a fifth-order Runge Kutta algorithm with adaptive step size 
[22]. 

With the grid algorithm described above the flow equations are solved 
starting the evolution in the broken phase deep in the UV region at a given 
k — A where we use a tree-level parameterization of the potential with only 
two initial parameters. For the results presented below the values A = 800 
MeV, = (140 MeVy, Aa = 49 and Za = 0.8 have been chosen. 

The evolution with respect to the scale k of the dimensionful minimum 
of the potential $o,fc in units of ^o(T = 0) and the quartic coupling for 

= 4 is shown in Fig. 1 for different temperatures without wavefunction 
renormalization corrections^''. The initial values at the UV scale A are fixed 
at T = and are kept constant for all temperatures. This, in principle, 
limits the predictive power of the high-temperature behavior of the calculated 
quantities. We have checked that for the range of temperatures relevant for 

^^Including wavefunction corrections does not change the shown curves significantly. 
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In (k/Oo) In (k/$o) 

Figure 1: The fc-evolution of the minimum of the full potential 
(left panel) and of the quartic coupling (right panel) for various 
temperatures. 

the phase transition the results are insensitive to temperature variations in 
A within reasonable limits. 

For all temperatures below the minimum $o,fc/^o(2^ = 0) (left panel) 
decreases as a function of k and settles down to a constant for very small 
/c-values where we can stop the evolution to obtain the renormalized vacuum 
expectation value ($o) (cf- [23]). Which is the physical order parameter for 
the 0{N) symmetry phase transition. 

The quartic coupling which in four dimensions is a dimensionless quan- 
tity shows a logarithmic running and tends to zero for A; — as can be seen 
in the right panel of Fig. 1. For > 1 there are (A^ — 1) massless pions 
in the spontaneously broken phase. Below the phase transition they never 
decouple from the evolution and thus always contribute to the logarithmic 
running resulting in an infrared-free theory. The undulation in the curve for 
small temperature (T = 10 MeV) is a numerical artefact and stems from the 
finite grid spacing. 

4.2 The potential and the quark condensate 

In the infrared the order parameter $o signals a second-order phase transition 
(see Fig. 2). The vacuum expectation value is obtained as the dimensionful 
minimum of the full effective potential with respect to the field, including 
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Temperature [MeV] 



Figure 2: The order parameter for the 0(4)-model as a function of 
temperature. A comparison with the chiral perturbation theory ex- 
pansion for low temperatures and with the Widom scahng function 
in the vicinity of the critical temperature is also shown. 

wavefunction renormalization contributions. For low temperatures we com- 
pare our result for = 4 with a three-loop chiral perturbation theory expan- 
sion up to order O {(T/ f.„)^) [24] and find perfect agreement up to T !^ 45 
MeV. The critical temperature is a non-universal quantity and therefore de- 
pends on the initial UV parameters. Our value of ~ 140 MeV is higher 
than the value ~ 100 MeV of ref. [6] also obtained for a full potential 
analysis with the exact renormalization group including wavefunction renor- 
malization contributions. This discrepancy in the value of Tc might be due 
to the different implementations of the threshold functions in the flow equa- 
tions. In the vicinity of we see a scaling behavior of the order parameter 
with a critical exponent f3 ~ 0.4 also depicted in Fig. 2. 

In our analysis we are not restricted to universal quantities such as the 
minimum of the potential. Indeed we do know the entire potential as function 
of the field $ which is shown in Fig. 3 (left panel) in the — > limit for 
various temperatures [T = 0, T = and T > Tc). Here one again observes 
a second-order phase transition. Numerically we have to stop the evolution 
at a very small but finite fc-value which slightly influences the final shape of 
the potential. In fact, the potential at = corresponds to the free-energy 
density which should be a fiat function around the origin due to convexity. 
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For a finite A;- value we always obtain the well-known Mexican hat form for the 
broken phase (right panel). The evolution towards the Maxwell construction 
([25]) at T = can explicitly be seen. 




Figure 3: The dimensionful convex potential for different tempera- 
tures (left panel) and the fc-evolution of the dimensionful potential 
at T = (right panel) towards zero. 



5 Critical behavior 

Within the RG approach we can directly link the zero-temperature physics 
to the universal behavior in the vicinity of the critical temperature. We 
will restrict ourselves to temperatures T < Tc in this section but there is in 
general no restriction to the broken phase. 

In order to investigate the critical regime of the system it turns out that 
the use of rescaled quantities is convenient since the scale k is no longer 
present explicitly in the evolution equations. At the critical temperature only 
properly rescaled quantities asymptotically exhibit scaling and the evolution 
of the system is purely three-dimensional. This is the well-known dimensional 
reduction phenomenon (see e. g. ref. [26]) and is exhibited in the structure 
of the above mentioned threshold functions. It is also the reason for the 
fractional power in the threshold functions at finite temperature in Eqs. (22)- 
(24). For the dimensionally reduced system at high temperature (T ^ k) 
only the bosonic zero-mode Matsubara frequencies are important while all 
the other modes decouple completely from the evolution. 
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We define the rescaled (dimensionless) quantities as follows 

02 .= j^-{d-^+v)^2 ^27) 



where we have introduced the anomalous dimension 

?7 = -^lnZfc. (28) 

As already mentioned in Sect. 3 there are in general two anomalous dimen- 
sions for the 0(A^)-model, one for the radial mode and another one for the 
Goldstonc bosons. Due to our uniform wavefunction renormalization approx- 
imation wc work with only one rj^^. 

In this section all results are obtained by taking the corrections due to 
the anomalous dimension rj explicitly into account. In order to deal with 
dimensionless flow equations we additionally rescale the potential via 

u{(j)^) := k-''V{(j)^) . (29) 

All dimensionless quantities are denoted in the following by small letters. 

From Eqs. (19)-(20) we obtain the following properly rescaled dimension- 
less flow equations in d dimensions 

c 

dtu{(j)'^) = -du + {d - 2 + 7])(i)\' + ^ 



1 



+ 



A^- 1 



l + 2u' + A(f)^u" 1 + 2u' 



(30) 



2S, 



1 + 



+ 



(1 + 402ti")2 202^" \l+ A(t)^u" 



<t>=<t>0 



where the primes on the potential u indicate differentiation with respect to 
(fp'. Since no scale k appears explicitly in the rescaled equations they are in a 
scale-independent form. Due to the dimensional reduction phenomenon we 

do not need to use the finite temperature flow equations at the critical tem- 
perature. It is sufficient to employ the three-dimensional zero-temperature 
equations in oder to investigate the critical regime of the phase transition. 



'^At the scaling solution both anomalous dimensions become degenerate anyway. 
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5.1 The scaling solution and fixed points 

We solve the set of coupled equations (30), (31) for arbitrary N with the 
same initial values at t = 0^^. We arbitrarily set A = 0.5 which is anyway 
an irrelevant quantity and finetune only 0olt=o at t = in order to find a 
fc-independent solution, the so-called 'scaling solution' of the flow equations. 
A second-order phase transition involves an infrared fixed point of the renor- 
malization group. Therefore the physics close to the phase transition is scale 
invariant and the critical behavior should be described by a scale- independent 
solution. 

For an initial value around the critical 0oltlo scale t = this 

is indeed the case and the evolution towards the infrared t — — oo reaches 
almost the scaling solution where the t-dependence vanishes and 0q tends to 
a constant (fixed point) value (f)^* . This feature is demonstrated in Fig. 4 
where the evolution of (pQ versus the "flow time" t = ln(A;/A) is shown for 
different initial values at the UV scale t = 0. 




ln(k/A) 



Figure 4: The evolution of the dimensionless minimum 0q as func- 
tion of the "flow time" t = ln(A;/A) for = 4 and different initial 
values at t = around the critical value 0oltlo (^^^ ^^-^^ ^"-"^ details). 

For a starting value near the critical value (plltLo the evolution deviates 
either towards the spontaneously broken (0q ^ 0) or the symmetric (0q = 
0) phase in the infrared as shown in Fig. 4. Due to the rescaling of the 

^^For the rescaled system no explicit value for the UV cutoff A is necessary. 
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dimensionful minimum $q (cf Eq. (27)) which in the infrared tends to a 
constant value or zero during the evohition with respect to the scale k, the 
dimensionless minimum 0q = $o/^^ diverges for k ^ for the broken phase 
as can be seen in Fig. 4. 

In the further course of the evolution, for a initial value very close to 
the critical one at t = 0, the scaling regime is approached where the system 
spends a long "time" before it finally deviates from it again. The scaling 
solution for the minimum of the potential is denoted by a straight line in 
Fig. 4. 

Not only 0q but all other quantities reach the scaling solution around 
t !v —10 and leave the critical trajectory at t —18. The time that the 
system spends on this scaling solution can be rendered arbitrarily long by 
appropriate finetuning of the initial values at the UV cutoff scale i.e. at t = 0. 
In order to produce such a high precision as shown in the Figure one has to 
know the initial starting value up to 13 digits. 

Besides of the trivial high-temperature and low-temperature fixed points 
the 0(iV)-symmetric model exhibits a nontrivial mixed fixed point inherent 
in the fiow equations. The region around this fixed point is mapped out 
in Fig. 5 where the rescaled dimensionless quantities (pl and A are plotted 
for = 4 and for different initial values during the evolution towards the 
infrared. The arrows in the Figure characterize the flow w.r.t. the scale t 
towards zero. 

In general, there are two relevant physical parameters for the 0{N) model 
which must be adjusted to bring the system to the critical fixed point accord- 
ing to universality class arguments. Consider for example the ferromagnetic 
spin Ising model with a discrete symmetry Z2 — 0{1) where the relevant pa- 
rameters are the temperature and external magnetic field. Since we work in 
this section in the chiral limit (no external sources or masses) only one rele- 
vant eigenvalue from the linearized RG equations is left, the temperature^"^. 
This can be seen in Fig. 5. The quantity (pl is the relevant variable because 
repeated renormalization group iterations (which are the discrete analog to 
the continuous evolution with respect to t — > —00) drive the variable away 
from the fixed-point value while A is the irrelevant variable and is iterated 
towards the fixed point if the initial values are chosen sufficiently close to 
the fixed point. Thus one has a one-dimensional curve of points attracted to 
the fixed point, the so-called 'critical surface' where the correlation length 

^^The temperature-like relevant variable is also called thermal scaling variable. 
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diverges. 

0.18 
0.16 
00.14 

-e- 

0.12 
0.1 

1.8 2 2.2 2.4 2.6 2.8 

I 

Figure 5: The mixed fixed point for = 4. 

One nicely sees that the critical surface (line) separates both phases. 
When we choose initial values near the critical line the systems spends a long 
time near the critical point. For example starting with 0olt=o ^ ^o\t=o 
X finite, the evolution tends ultimately to the zero-temperature fixed point 
corresponding to very large values of 0q. On the other hand, starting below 
the critical line the evolution of 0g tends towards the high-temperature fixed 
point corresponding to small values of 0o- Only for exactly 0oL=o ~ ^o\t=o 
the renormalization group trajectory fiows into the critical mixed fixed point, 
which means that the large distance behavior at the critical point is the same 
as that of the fixed point. 

Here the predictive power of the nonperturbative approach becomes visi- 
ble: During the evolution near the scaling solution the system loses memory 
of the initial starting value in the UV and the effective three-dimensional dy- 
namics near the transition is completely determined by the fixed point and 
hence independent of the details of the considered microscopic interaction at 
short distances. 

The equation of state drives the potential away from the critical tem- 
perature. As a result after the potential has evolved away from the scaling 
solution its shape is independent of the choice of the initial values such as 
e.g. A for the classical theory at the UV, signaling universal behavior near 
the critical point. 
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5.2 Beta functions 

In Fig. 6 the beta function Px = dtX is presented as a function of A. For 
= 4 it has two zeros at A* = and A* = 2.28 corresponding to two fixed 
points. The position of the fixed points depends on the particular form of 
the RG fiow equation encoded in the choice of the blocking function fl^\ 
Physical results, however, such as critical exponents must not depend on a 
particular choice. 




-1 I ■ ' ■ ' ■ 1 

12 3 

X 

Figure 6: The beta function f3x{X) for the 0(4)-model. The RG 
result is compared to an one- and three-loop e-expansion for three 
dimensions. The arrows denote the fiow w.r.t. the scale k towards 
the infrared. 

The anomalous dimension t] is nonzero at the nontrivial fixed point as 
shown in Fig. 7. In general it is difficult to calculate the zeros of the beta 
function since this requires knowledge of the physics beyond perturbation 
theory. The ultraviolet-stable but infrared-unstable fixed point at A* = is 
called the Gaussian fixed point because the fixed-point Hamiltonian yields 
a (simple) Gaussian partition function. The properties of this fixed point 
depend on whether the dimension d is greater or less than 4. When d is 
(slightly) smaller than 4 A is (slightly) relevant at the Gaussian fixed point 
and thus can fiow to another nearby fixed point at A* ~ 2.28. This point 
is called the Wilson-Fisher fixed point and the fiow with respect to k is 
indicated in Fig. 6 by arrows. Thus in the limit k —>■ the coupling A^ is 



21 



driven to this infrared-stable fixed point because the derivative of the beta 
function is positive. 




Figure 7: The anomalous dimension r/ as a function of the quartic 
coupling A. The flat line indicates t] at the physical vacuum. 

We can also verify the large- limit: The fixed point A* tends to zero for 
iV — >■ oo merging with the Gaussian fixed point (see Table 1 and Fig. 8) [27]. 
When the two fixed points are sufficiently close to each other it is then 
possible to deduce universal properties at one fixed point in terms of those 
at the other one. 

Universal quantities can be calculated by perturbation theory at fixed 
dimension or by the e-expansion where the parameter e = 4 — c? is related 
to the number of dimension of the system. The computation of the criti- 
cal exponents via the e-expansion beyond order 0{e) requires knowledge of 
complicated higher-order terms in the perturbative renormalization group 
equations. It is believed that the e-expansion (including the assumption of 
Borel summability) is a good approximation scheme if a sufficient number of 
higher terms are taken into account. At the end of the calculation one has 
to set e equal to one (or two) to obtain results of direct physical significance. 
One should bear in mind, however, that for a given case the e-expansion might 
be an asymptotic expansion giving eventually inaccurate results if truncated 
at too high an order. 

In Fig. 6 we compare the beta function j3x for N = A obtained numerically 
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Figure 8: The beta function /3a (A) for 
various values of demonstrating the 
large- behavior. 





0r 


A* 


N= 


=1 


0.062 


3.20 


N= 


=2 


0.087 


2.88 


N= 


=3 


0.115 


2.58 


N= 


-A 


0.144 


2.29 


N= 


=10 


0.336 


1.26 


N= 


=100 


3.369 


0.15 



Table 1: In three-dimensions 
the Wilson-Fisher fixed point 
tends towards zero in the large- 
iV hmit. 



within our nonperturbative RG approach with e-expansions for different loop 
orders. 

The e-expansion for an 0(A^)-symmetry up to three- loop order [28] is 
given by (for e = 1) 



/3a 



(AT + 8) 2 (3Ar+14) 3 

-q ^ -Q - -Q 

y 6 12 ^ 



33A^2 ^ 922Ar + 2960 + 96(5iV + 22)C(3) 
123 



with (^(3) = 1.202057 and the short-hand notation 

12A 



9 



(4vr)fr(|) 



At three-loop order we find agreement with our result up to A ~ 0.55 
while the one- loop order deviates from our result around A ~ 0.45. Of 
course, the e-expansion, despite of a failure from a numerical point of view, 
provides an unambiguous classification of the fixed point consistent with our 
nonperturbative RG results. 
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5.3 Critical exponents and scaling relations 

In the vicinity of the critical point one obtains a scahng behavior of the system 
which is governed by critical exponents. They parameterize the singular 
behavior of the free energy near the phase transition. Altogether there are 
six critical exponents a, (3, 7, 6, v and ri for the 0(iV)-model but only two 
of them are independent due to four scaling relations: 

a — 2 — du 

7 = (2-r/)z^ (32) 
^ _ d + 2-r) 
~ d-2 + ri ' 

In general, these universal critical exponents depend only on the dimension- 
ality of the system d and its internal symmetry. The exponents rj and S 
describe the behavior of the system exactly at the critical temperature and 
the remaining ones parameterize the system in the region of the fixed point 
around Tc (T 7^ Tc). 

In this work we have calculated the critical exponents f], f3, v and 8 in 
three-dimensions numerically for different values of N and could verify the 
scaling relations (32) with deviations typically being less than 0.1 %. The 
smallness of the anomalous dimension r] (see Table 2) justifies the use of the 
derivative expansion of the effective action. 

The critical exponent parameterizes the behavior of the spontaneous 
magnetization or the order parameter for the broken phase (0olt=o ^oltlo) 
in the vicinity of the critical temperature (see Fig. 2). It is here numerically 
obtained as the slope of the logarithm of 0q as a function of ln(0olt=o ~ 
0oltlo) ^'-'^ ^^^y small arguments. The difference (0olt=o~ ^oltlo) ^ measure 
of the distance from the phase transition irrespective of any given value 
At=o at the UV scale. If 0olt=o interpreted as a function of temperature 
this difference is proportional to the deviation from the critical temperature 
i.e. proportional to (T — T^. Thus the quantity 0oltlo defines the critical 
temperature in three-dimensions. 

The critical exponent v describes the scaling behavior in the critical region 
of the renormalized mass which is the inverse of the correlation length ^ = 
{T — Tc)~^ and can be defined for arbitrary N by the relation 



^l\t=o - ^l\T=o (33) 
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Critical Exponent v 



v = 0.645 



In \(j)l\t=o 



L2\cr I 



Figure 9: The critical exponent u for = 4. 



with the renormalized mass in the symmetric phase = k'^u'{0) for k 0. 
Wavefunction renormalization corrections are included in the potential. Ex- 
actly at the critical temperature i.e. for {(f>l\^^Q = ^oltlo) renormahzed 
mass vanishes and hence the correlation length diverges. For the determina- 
tion of the exponent u we calculate for different values of 0olt=o aground 



^o\t=o) 



^o\t=o 



The correlation 
Note that the 



and plot In(m^) as a function of In 
is linear (see Fig. 9) and the slope yields the exponent u 
definition of the renormalized mass results in a positive value only in the 
symmetric phase. Therefore one has to make sure to reach the symmetric 
phase via 4>lit) — > before the evolution with respect to t is stopped. 

In statistical physics the exponent 6 describes the variation of the mag- 
netization (here $) in an external field B in the limit B ^ at T^. 



B 



dV{<l>) 
9$ 



(34) 



This means that the rescaled potential ^(0^) is proportional to 0"^"*"^ for small 
fields. In mean-field theory where no quantum fluctuations are taken into 
account the potential near the origin starts with a 0^-term in four-dimensions 
at the critical temperature. This behavior changes completely in three di- 
mensions. In order to determine the critical exponent S we study the flow 
equations in the limit 0^ — > oo corresponding to k —>■ for fixed (not rescaled) 
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77 


5 


V 


/5 


RG(1) 
KG (2) 
N=0 lattice 
e 

SPS 


0.039 
0.034 
0.041(25) 
0.032(3) 
0.027(4) 


4.77 
4.80 
4.77(14) 
4.81(2) 

4.84(3) 


0.59 

0.592(3) 
0.588(2) 

0.588(2) 


0.30 

0.308(6) 
0.304(2) 
0.302(2) 


RG(1) 
RG(2) 
N=l lattice 
e 

SPS 


0.044 
0.045 
0.044(31) 
0.038(3) 
0.032(3) 


4.75 
4.74 
4.75(17) 
4.78(2) 
4.81(2) 


0.64 

0.631(2) 
0.631(2) 
0.630(2) 


0.32 

0.329(9) 
0.327(2) 
0.325(2) 


RG(1) 
RG(2) 
N=2 lattice 
e 

SPS 


0.043 
0.042 
0.038(18) 
0.040(3) 
0.033(4) 


4.75 
4.76 
4.78(10) 
4.77(2) 
4.81(3) 


0.69 

0.674(2) 
0.671(5) 
0.670(2) 


0.35 

0.350(5) 
0.349(4) 
0.346(2) 


RG(1) 
RG(2) 
N=3 lattice 
e 

SPS 


0.041 
0.038 
0.038(17) 
0.040(3) 
0.033(4) 


4.76 
4.78 
4.78(10) 
4.77(2) 
4.81(3) 


0.74 

0.711(2) 
0.710(7) 
0.705(3) 


0.37 

0.369(5) 
0.368(4) 
0.365(3) 


RG(1) 

N=4 ^^^^^ 
lattice 


0.037 

n no /I 
U.Uo4 

0.0254(38) 

osn 


4.79 
4.0U 

4.851(22) 


0.78 

0.7479(90) 
73('9'l 


0.40 

0.3836(46) 


N-7 ^^(1) 

' RG(2) 


0.027 

0.023 


4.84 
4.87 


0.86 


0.43 


RG(1) 
^ RG(2) 


0.021 
0.017 


4.88 
4.90 


0.91 


0.45 


N-lOO ^^^^^ 


0.0025 
0.002 


4.99 
4.99 


0.99 


0.49 


large- 





5 


1.0 


0.5 



Table 2: Critical exponents for different iV compared with a perturbative 
calculation (SPS), the e-expansion (e) [28] and lattice results [29]. RG(1) 
and RG(2) denote our results. See text for details. 
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In this limit u' also increases and we can therefore simplify the flow 
equations by neglecting the threshold functions (see Eq. (30)). At the scal- 
ing solution dtu = (corresponding to T = Tc) we can solve the equation 
analytically to find 

yielding the well-known scaling relation between S and r) (see Eq. (32) and 
[23]). 

In order to extract the exponent 6 numerically we calculate the first 
derivative of the potential u w.r.t. the rescaled field (f)^ and determine the 
slope of Inu' for large In^^. 

Table 2 summarizes the resulting critical exponents r], S, v and /3 for 
various iV within our RG approach (without quoting small numerical errors). 
The label RG(1) denotes the numerical results for the exponents r/, v and /5 
which we have obtained from the scaling solution (cf. e.g. Fig. 4). 

Via the scahng relations (32) we can, on the other hand, directly calculate 
b once knowing 77. The abbreviation RG(2) labels the calculated result for the 
exponent 77. We also compare these exponents with a lattice study [29], an 
expansion (denoted by e) and to a Borel-resummed perturbative calculation 
up to three-loops at fixed dimension three (labeled by SPS) [28]. 

One nicely observes the convergence of all calculated critical exponents 
to the large- A?^ values which are quoted in the last line of Table 2. It is known 
that in the large-A^ limit, the local potential approximation becomes exact 
and the anomalous dimension 77 vanishes [9] . For the direct determination of 
the remaining critical exponent a we would have to calculate the specific heat 
at the critical temperature which was done in [5] . We omit the quotation of 
this exponent in Table 2. 

5.4 The influence of the blocking functions 

The structure of the regularized RG-improved flow equations (but not the 
universal results) depends on the a priori unknown cutoff function jf'. The 
index i refers to the order of derivatives (9^*) in the gradient expansion of 
the action. Thus blocking functions with i > \ regularize higher-order con- 
tributions to the action which go beyond the uniform wavefunction renor- 
malization approximation, used here. 

In a previous paper [4] we have studied different choices of the cutoff 
function for a truncated set of flow equations (A similar analysis was done in 
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[30]). In order to see this influence for the non-truncated set of flow equations 
for the full potential we derive flow equations for a similar choice of the cutoff 
functions as in ref. [4] and calculate some critical exponents. 

To be consistent with the structure of the higher momentum integrals we 
evaluate differential equations for each higher blocking function. Without 
the explicit solution of these blocking-function equations we can still find the 
corresponding set of flow equations (without the wavefunction contribution). 
We list below the corresponding flow equations for cutoff functions i = 0, 1, 2 
although fj^^ is not needed in our approximation scheme: 



Ik 



(1) 



(2) 



OtV = 



+ 



N-1 



(l + m^/F) (l + m2/F) 

1 N-1 



d{d + 2) 



(l + m2/A;2)2 {l + ml/k^ 



(36) 
(37) 

(38) 



As an example we take N — 4 and compare the dependence of the critical 
exponents u and (3 on the blocking function of order i. 



i 


V 


/5 





0.853 


0.42 


1 


0.815 


0.405 


2 


0.814 


0.403 



Table 3: The critical exponents v and /3 
for N = 4 with different blocking func- 
tions Z^*'' (without wavefunction renor- 
malization) . 



In Table 3 the results are listed for various cutoff functions . A small 
systematic decrease in the values with the order i is observed. Already for 
i > 1 there is almost no difference which validates the findings in [4]. This 
permits the conclusion that for higher blocking functions the convergence 

^^See notation and definitions in Section 3.2. 
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becomes more and more stable and a decrease towards the values found by 
lattice simulations and e-expansions is observed. 

6 Summary and conclusions 

Near a phase-transition a perturbative treatment fails as an adequate de- 
scription of the critical behavior and nonperturbative methods such as the 
RG method become necessary. We have presented a RG study for the finite- 
temperature phase transition of self- interacting scalar theories with 0{N)- 
symmetry. We have used an implementation of the Wilsonian approach to 
field theories in thermal equilibrium which provides a direct link of the 0{N) 
universal behavior near the critical temperature with the physics at zero 
temperature. It offers a nonperturbative way of studying the mechanisms of 
dimensional reduction. 

Based on an application of Schwinger's proper-time (heat kernel) regu- 
larization within a one-loop expansion of the effective Euclidean action we 
have introduced various cutoff functions by comparing the proper-time reg- 
ularization with a sharp-momentum regularization rendering the resulting 
flow equations infrared finite. This specific choice of cutoff functions has the 
additional advantage that all resulting flow equations can be expressed in a 
transparent and analytical form. 

We have generalized previous studies [5], [4] to a general 0(A^)-symmetry 
by going beyond the local potential approximation in a derivative expansion. 
We have solved the highly coupled set of flow equations numerically on a grid 
without an specific ansatz or truncation for the potential and its derivatives 
during the evolution. In the infrared a convex potential is obtained for all 
temperatures, consistent with the concept of the effective action approach 
which is based on a Legendre transformation. 

At the critical temperature all dimcnsionful quantities can be rescaled 
which simplifies the numerics. The dimensional reduction phenomenon which 
is embedded in a very transparent way in our threshold functions could be 
confirmed. We have discussed the critical behavior of the 0{N) model in 
detail, by determining fixed points and beta functions for various N. We 
have calculated independently four critical exponents verifying the scaling 
relations among them. 

One important ingredient in our approach is the choice of the blocking 
function which is a priori unknown and could infiuence the universal quanti- 
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ties via the convergence properties of the flow equations. In order to control 
this influence we have recalculated two critical exponents (z/ and f3) for = 4 
for three different choices of the blocking functions omitting the wavefunction 
renormalization. No strong dependence has been observed. Within numer- 
ical errors we find stable exponents except for the lowest blocking function 
fj^^^ due to a poor convergence. For higher blocking functions the convergence 
becomes more and more stable in agreement with [4]. 

Altogether the presented RG approach has proven to be a powerful non- 
perturbative tool in the study of finite-temperature field theory. Due to the 
omission of any temperature dependence of the initial conditions in the ultra- 
violet its reliability should be best for low temperature. For N — A our results 
indeed agree perfectly with chiral perturbation theory. In contrast to chiral 
perturbation theory we can, however, extend the temperature range to the 
critical temperature where the phase transition takes place. This encourages 
the application of this approach to other theories involving e.g. gauge-fields. 
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A Appendix 

In order to improve the readability of Sec. 3 we collect in this appendix some 
technical details concerning the derivation of the wavefunction renormaliza- 
tion flow equation (12). 

The expansion of the effective 0(iV)-invariant Lagrangian to two-deriva- 
tive order (see Eq. (7)) requires the calculation of a trace over (N x iV)-matrix 
valued terms stemming from different powers of the second derivative of the 
potential V/j with i,j = l,...,N. 

This can basically be accomplished by introducing projectors 



and X2 = Sij — 



$2 



$2 
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which satisfy the usual projection operator property XaXj, — 5abXi, for a,b — 
1, 2 and trXi = 1, trXa = iV - 1. 

Using the abbreviations ml = A(3<l^ — $0) ^"^^ ^1 = ~ ^0) 

(cf. Eq. (8)) the N x N matrix V/j can be rewritten by 

y/' = A - ^l) 5ij + 2X^i^j = mlX^ + mlX, . 

Using the projector properties different powers of V-j can be readily evalu- 
ated. One finds, for example, 



Tr 



E 

Ln=2 



1 \n—l^n n—1 



k=0 



'-Tr [X2 



(A.l) 



^ml - ml 



Tr 



x.dlx. 



where on the l.h.s. of Eq. (A.l) all finite and infinite sums are analytically 
calculable (cf. [17]) and the remaining traces are given by a straightforward 
calculation as 



Tr[Xad^X, 



Tr 







with a,b,c = 1,2 . 



XadlX,] = 2(-l)'^<"'+l 



$2 



Tr[Xad^X,d^X,] = (-1)'^ 



V 



$2 



Evaluating the other contributions to the effective action in the same manner 
yields Eq. (10). 
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